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Abstract. Bose-Einstein correlations in relativistic heavy ion collisions and 
their dependence on the freeze-out condition in hydrodynamic models is ex- 
amined. The Cooper & Frye space-like freeze-out mechanism is compared to 
time-like freeze-out, where particles are emitted away only from the surface, 
i.e. space- vs. time-like freeze-out. The corresponding HBT radii are cal- 
culated for the two models emphasizing the difference in the outward HBT 
radius. 



1. Introduction 

Bose-Einstein interference of identical particles or the Hanbury-Brown & Twiss ef- 
fect (HBT) jl| shows up in correlation functions of pions and kaons emitted from the 
collision zone in relativistic heavy ion collisions. It is an important tool for determin- 
ing the source at freeze-out and recent data from relativistic heavy ion collisions can 
restrict the rather different models, that have been developed to describe particle 
emission in high energy nuclear collisions. In hydrodynamical calculations particles 
freeze-out at a hypersurface that generally does not move very much transversally 
until the very end of the freeze-out HH]. In cascade codes the last interaction 
points are also found to be distributed in transverse direction around a mean value 
that does not change much with time [^-^], but the width of the emission zone 
increases from narrow surface emission to a widespread volume emission. In the 
first stage the freeze-out surface is relatively static in the transverse direction and 
the emission extends in time whereas in the second stage the freeze-out happens 
relatively rapid all over the spatial extent of the source. The two stages are often 
referred to as surface and volume freeze-out respectively, or time-like and space-like 
freeze-out (see Fig. l). a The distinction is, of course, only approximate and the 
amount of particles assigned to the two stages varies between models. 

In the Cooper & Frye freeze-out assumption the net number of particles leaving 
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Fig. 1. Schematic picture of the transverse position and invariant time of the 
freeze-out surface in hydrodynamic calculations. The dashed curve approximates 
the source by time- and space-like components. 



the hydrodynamic region is fic| 

Cooper & Frye (space — like) : = — [ da ^p^ f (x , p) : (1) 

d A p 2tt J a 

where f(x,p) — (exp(p- u(x)/T(x) ± 1) _1 is the standard thermal Bose or Fermion 
distribution function with local flow velocity u(x). The freeze-out hypersurface is 
denoted by and is usually defined at a constant energy density or temperature 
T c ~ 100 — 150 MeV which is determined from the p± slopes. 

It was pointed by Sinyukov |ll| that the Cooper & Frye freeze-out was inconsis- 
tent for time-like freeze-out since particles might reenter the hydrodynamic phase. 
He therefore reflected the momenta in the local rest frame at the boundary for those 
particles that would have reentered the fluid 

Sinyukov (time - like) : E^- = A J da^f{x, K a [p]) , (2) 

where lZ a [p] reflects the component of the particle momentum normal to the freeze- 
out surface in the local rest frame of the fluid. The argument for doing this was that 
for each particle reentering there is an outgoing particle inside the fluid that would 
have escaped into the free streaming region. However, this is not true when the 
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surface is moving (with velocity v s ) in which case slow particles may be overtaken 
by the freeze-out surface. Defining n CT (x) as a vector normal to the freeze-out surface 
and directed away from the fluid and outward into the free streaming phase, the 
condition for freeze-out is that particle velocities v p = p/E obey 

(v p -v s (x))-n ff (x) >0. (3) 

Particles in the fluid cell to be frozen out that do not fulfill this condition should 
therefore have their momenta reflected in the frame of the surface. When the surface 
moves rapidly inwards not many particles need to have their momenta reflected. 
Likewise, when there is strong outward directed flow the thermal factor exp(p-u(x)) 
in distribution function automatically guarantees that not many particles move 
slower than the surface speed. But for small flow and slowly moving surfaces half 
of the particles should be reflected. 

The above prescription for correcting the time-like part of the Cooper & Frye 
freeze-out assumption may be considered as an improvement. Yet, there are a 
number of considerations that remain uncured. It should be pointed out that this 
modified time-like freeze-out assumption preserves momentum and energy globally 
in the fluid, when it is applied to both sides, but not locally near the surface. At 
the surface there is a net back reaction due to the surface emission/evaporation 
that in the time-like freeze-out is away from the surface only. In the treatment 
of Godunov (see, e.g., [|l2|]) restricting freeze-out to those particle for which the 
four-vector product, n(x) ■ p > 0, is positive does, however, not conserve global 
energy. Also rescatterings and final state interactions of emitted particles are still 
ignored. Another issue is the assumption of a sharp freeze-out surface. In reality 
the freeze-out take place when the particle mean free path becomes long enough 
that it can escape and the hydrodynamic assumption of a short mean free path 
breaks down near the surface. This is naturally included in cascade models which 
find an extended region of final interactions points around what is thus a diffuse 
surface layer. Recently, Grassi et al. have improved the sharp surface assumption 
in hydrodynamic models by a Glauber model approach ]l3] ] (see also Csernai, these 
proceedings Q). Similar ideas have been discussed in [|l4j in connection with opaque 
sources emitting from a surface layer of thickness ~ A m / P and their effects on HBT 
have been calculated. 

2. Space-like vs Time-like emission 

We will discuss space-time correlations from hydrodynamic sources with the Cooper 
& Frye (space-like) and Sinyukov (time-like) freeze-out assumption and compare 
them. In the subsequent chapter we will then consider two-component sources with 
both space- and time-like emission. 

To study the differences between the freeze-out conditions we restrict ourselves 
to a class of models that have: 

i) cylindrical symmetry around the beam or z-axis, b 
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ii) longitudinal expansion with Bjorken scaling, 
Hi) no transverse flow. c 

The distribution of emitted particles also referred to as the source of emission 
points is with the Cooper & Frye freeze-out assumption 



Here, the Bjorken variable r = \/t 2 — z 2 is the invariant time and r\ — 0.5 ln(i + 
z) / (t— z) the space-time rapidity. The flow four- vector is u = (cosh(7/), sinh(r/), 0, 0), 
which gives p ■ u = m± cosh(?7 — Y). R(t) is the transverse radius of the freeze-out 
surface that moves in time. Initially, when the nuclei collide R(0) is the transverse 
size of the nuclear overlap zone and final freeze-out takes place when R(rf) = 0. 
The surface speed is v s = dR(r)/d,T. The temporal source factor S t (t) determines 
the amount of particles that freeze-out per surface element at invariant time r. 
Notice that any normalization is irrelevant as they cancel out in later correlation 
functions ( |l5|) and HBT radii. 

The Sinyukov freeze-out assumption for a time-like source without flow differ 
from the Cooper & Frye by 



where 9 is the polar angle between p and the vector n a (x) normal to the freeze-out 
surface and directed away from the fluid and outward into the free streaming phase. 
Notice, that the reflection does not change the single particle distribution functions 
but only correlation functions that are sensitive to correlations between p and x. 



We follow the standard definition of correlation functions and the HBT radii as is 
described in more detail in the appendix. It is common to boost longitudinally so 
that the pair rapidity always is Y = and choose the direction of p = (pi + P2)/2 
along the outward or x direction whereas the third y-direction is called the sideward 
directions. In cylindrical coordinates we choose the polar angle with respect to the 
x-axis and so the factor distinguishing the Cooper & Frye from the Sinyukov freeze- 
out is simply 9(cos#). 

The HBT radii are for a space-like source Scf(p, x), for which (x) = (cos 9) = 0, 



S CF (x,p) ~ e-^*)/n*) S T (r)S(r ± - R(t)) , 



(4) 



S s (x,p) = 2S CF {x,p)e{p ■ n CT (x)) = 2S C f{x,p) 9(cos0) 



(5) 



3. HBT Radii 






<(z-Ai) 2 )^(r 2 )— . 

(y 2 ) = \(R(Tf) , 

((x-/3 t) 2 ) = i(i?(r) 2 ) +P 2 a(r). 



(8) 



(7) 



(G) 



The temporal and angular averages and fluctuations are described in the appendix. 
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The outward HBT radius is thus larger than the sideward [|l6|,|l7|] 

R 2 o = R 2 s +0 2 o a(r). (9) 

The measured out- and sideward HBT radii are very similar in relativistic heavy 
ion collisions and is based on Eq. (|^) taken as an indication of very small duration 
of emission, i.e., particles appear or freeze-out in a "flash" Jl9[ . 

Changing the Cooper & Frye freeze-out to that of Sinyukov does not change 
the sideward and longitudinal HBT radii since the reflection is only in the outward 
or x-direction. However, the outward HBT radius differ because the reflection has 
the effect that only particles from the front side of the source pointing towards a 
given detector are measured in that detector whereas those from the back side are 
not. The geometry is chosen such that the outward or x-axis is always along the 
particle momenta p ~ p which guaranties that only the half sphere toward any 
detector is seen in that detector. This has the important consequence that (x) no 
longer vanishes but almost cancels (x 2 ) 1 / 2 . We fmd d 

Rl = {\-- 2 ){R{T)f+Pl<?{T) + \a(R(T)) 

I 7T Z 2 

-2/3 -((i?(r)-(i?(r)))(r-(r))). (fO) 

For a static surface R(j) — constant the last two terms vanish; for an inward moving 
surface the last term is positive. The first term in (QOj) is <r(x) and is significantly 
reduced as compared to (||) . As a consequence Eq. (J^j no longer holds as was also 
found for opaque sources or sources with strong transverse flow fl4| . 



4. Two Component Source 

It is important to distinguish between space- and time-like freeze-out, since they 
give very different HBT radii. Hydrodynamical models assume that parti- 

cles are emitted from the surface. This is actually also found in some cascade models 
at early times of the collision ^-^] , but eventually the whole source freezes out and 
disintegrates. The late stage of cascade models resembles more a volume freeze-out. 
These sources can approximately be described by two components, initially surface 
emission but eventually volume freeze-out. Generally, for a two-component source 
S(x) = pSi(x) + (1 — p)S 2 {x), properly normalized (J d 4 xSi(x) = 1) such that p is 
the fraction of particles from source 1, the fluctuations in a quantity O is from ( |l7| ) 

a(0) =p*i(0) + (1 - p) a 2 (0) +p(l-p) ({0) 1 - {0) 2 f . (II) 

Here, (O)i = J d 4 xSi(x)0 and a % {0) = (0 2 ). t - (O) 2 . The fluctuations are the 
weighted sum of the fluctuations in the individual sources and an additional cross 
term. Since (y) = and (z — (3it) also vanishes for Y = this additional cross term 
does not contribute to the sideward and longitudinal HBT radii. It is, however, 
nonnegligible for the outward HBT radius. 
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Let us approximate the qualitative features of hydrodynamic and cascade mod- 
els by a toy model with two-components: an approximately static surface R(t) = R 
for times < t < r/ and a rapid freeze-out of the remaining source at t = tj. 
The first component violates Eq. (|^) for the inward moving particles and is thus 
time-like for which we should use the Sinyukov freeze-out assumption. The second 
component is space-like and we should apply the Cooper & Frye freeze-out assump- 
tion. We weight the particles frozen-out from the two component by the fractions p 
and (1 — p) respectively. The corresponding HBT radii are found by the expressions 
(r)j ©' ©' an d (0) by combining them according to ( |TT| ) 



(1 ?u T 

?2 _ !/i x D 2 



Rl = d-ftf-. (12) 



Ri = -(l+p)R\ (13) 



^ = U\£+p)-^)# + &fii, (14) 

In Fig. 2 we show the dependence of the out- and sideward HBT radii on 
the space- vs. time-like fraction p. The freeze-out assumption clearly plays an 
important role and the difference R 2 — R 2 S is a very sensitive quantity. The fraction 
of time-like emission p is a model dependent parameter and in reality there will 
be a gradual transition from time- to space-like emission. Also a finite width of 
the emission layer |Q will diminish the difference. Transverse flow has a similar 
effect of directed emission and leads to further reducing R a with respect to R s Q . 
However, strong transverse flow u± also has the effect of improving the condition 
of Eq. (|^) reducing range of particle velocities for which time-like emission occur. 



5. Summary 

The validity of Cooper & Frye space-time freeze-out assumption has been discussed. 
For sources with time-like emission it breaks down and one should rather apply a 
modified freeze-out assumption like that of Sinyukov which only allows emission 
away from the fluid. It was shown in a simple toy model that the different freeze-out 
conditions gave drastically different results for the outward HBT radius. A number 
of effects such as transverse flow, diffuse surfaces, etc., reduce the difference between 
the two freeze-out conditions. 

In the hydrodynamic calculations of Ref. || at RHIC energies the outward 
HBT radius R a is considerably larger than R s . This is mainly due to the existence 
of a long lived mixed phase of quark-gluon plasma and hadronic matter for the 
equation of state and initial conditions assumed in this calculation. Thus it is the 
duration of emission that gives large R 2 — R 2 or R a /R s and this part is unaffected 
by changing the freeze-out assumption from Cooper & Frye to that of Sinyukov 
as can be seen from Eqs. (||) and ([l(]). At AGS and SPS energies, however, the 
measured R s and R s are very similar excluding the possibility for a long lived mixed 
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P 

Fig. 2. Sideward and outward radii as function of the fraction of time-like emitted 
particles p; the amount of space-like emitted particles by the standard Cooper & 
Frye freeze-out assumption is (1 — p). Dashed curves are Rq/R 2 with freeze-out 
times of (3 Tf/R = 0, 1, 2. 

phase and the difference arising from the freeze-out assumption may be significant. 



Appendix: Correlation functions and HBT radii. 

For the correlation function analysis of Bose-Einstein interference from a source 
of size i?, we consider two particles emitted a distance <~ R apart with relative 
momentum q = (pi — P2) and average momentum, p = (pi +p2)/2. Typical heavy 
ion sources in nuclear collisions are of size R ~ 5 fm, so that interference occurs 
predominantly when q^h/R ~ 40 MeV/c. Since typical particle momenta are pt ~ 
p ~ 300 MeV, the interfering particles travel almost parallel, i.e., p\ ~ P2 — p 3> q. 
The correlation function due to Bose-Einstein interference of identical particles from 
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an incoherent source is (see, e.g., O]) 



/ d 4 x S(x,p) e iqx 2 



a '*'> = 1± ' j**s(*.p) 1 ■ (15> 

where p) is a function describing the phase space density of the emitting source. 
The +/— refers to boson/fermions respectively. 

Experimentally the correlation functions for identical mesons (7r ± 7r ± , if ± if ± , 
etc.) are often parametrized by the gaussian form 

C 2 (q s ,q , qi) = l + Aexp^i? 2 - fin* - q fRf - 2q q l R 2 ol ) . (16) 

Here, q = ki — k2 = {q Sl q ,qi) is the relative momentum between the two parti- 
cles and Ri,i = s,o,l the corresponding sideward, outward and longitudinal HBT 
radii respectively. We will employ the standard geometry, where the longitudinal 
direction is along the beam axis and the outward direction is along p and the side- 
ward axis is perpendicular to these. Usually, each pair of mesons is lorentz boosted 
longitudinal to the system where their rapidity vanish, Y — 0. Their average mo- 
mentum p is then perpendicular to the beam axis and is chosen as the outward 
direction. In this system the pair velocity f3 p =p/E p points in the outward direc- 
tion with (3 a = pj_/m± where m± = U m 2 + p 2 ± is the transverse mass. As pointed 
out in jl7| the out-longitudinal coupling R a i vanishes to leading order when Y = 0. 
The reduction factor A in Eq. ( |l6| ) may be due to long lived resonances |16|, 
coherence effects, incorrect Gamov corrections or other effects. It is found to be 
A ~ 0.5 for pions and A ~ 0.9 for kaons. 

It is convenient to introduce the source average and fluctuation or variance of 
a quantity O defined by 

<°> s ' } %% P l) • " (0,s(0,| - (01 ' (17) 

With qx ~ q • x — q-/3 p t one can, by expanding to second order in qiRi and com- 
paring to Eq. (|l|), find the HBT radii r\, i=s,o,l. They are @ 

R 2 = a( Xl -pit). (18) 

The HBT radii are a measure for the fluctuations of (xi — Pit) over the source 
emission function S. 

In the local center of mass system defined by Y = we have (3i — and 
Po = P±/m±. The HBT radii reduce to (x = R(t) cos and y = R(t) sin 9) |0 

R? = ct(z) = cr(rsinhr/) ~ (r 2 )— . (19) 

Rl = a(y)^(R(rf)a(sm9) = ^{R(rf), (20) 

Rl = ct(x-/3 <) = (R(t)) 2 ct(cos 9) + a (R(t)) (cos 2 9) + p 2 a{r) 

-2p o (cos9)((R(r) - (R(r)))(r - (r))) . (21) 
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The averages simplify because the space-time rapidity, angular and temporal 
integrations separates and due to the normalization a number of factors cancel. For 
example, a function of proper time only needs to be averaged with respect to the 
temporal parts of the source 



ft drrR(T)S T (T)0(T ) 
J Q Tf drrR(T)S T (r) 



<0(r)> = J ° rf ^-^oZr (22) 



In the toy model employed above the emission per surface element S T (r) is assumed 
constant. 

The angular averages also simplify for cylindrical geometry. ^From the defini- 
tions in Eq. ( |l7| ) we obtain 

r /2 , 9 0(9) cos 9d9 
(O(0)) = ~^j- 2 ■ (23) 

f 12 cose dd 

J — 7T/Z 

Notice that always (y) = (sin 9) — 0, whereas (x) = (R(t)) (cos 6) ^ 0, when 
cylindrical symmetry is broken as for the Sinyukov freeze-out condition, opaque 
sources, or sources with transverse flow. With transverse flow u(rj_) the thermal 
factor leads to a factor exp(p • u(r±)/T) in Eq. ( p3| ) |1| which moves the measured 
emission region in the direction towards p or the detector. 
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Notes 



a. The terminology chosen here is to call surface emission time-like as it takes 
place from a small volume but for a long time, whereas volume emission is 
space-like as it occurs in a large volume in a short period of time. The opposite 
terminology is sometimes used when referring to the direction of the four- vector 
n r7 (x) normal to the freeze-out hypersurface. 

b. However, cylindrical symmetry is broken by the direction of the detector when, 
for example, the source is opaque, has transverse flow, or for time-like freeze-out. 

c. Transverse flow has been included in the case of opaque sources |l5[. We will 
comment on the effect of flow later. 

d. In case of time-like freeze-out, the factors 2/tt arise from the nonvanishing 
angular average (cos 9) = f*^, 2 cos 9 d6/ Jl^ 2 / 2 d9 = 2/w. 
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